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ABSTRACT 

Using a simple uniform jet model of prompt emissions of gamma-ray bursts 
(GRBs), we reproduce the observed peak energy-isotropic energy relation. A 
Monte Carlo simulation shows that the low isotropic energy part of the relation 
is dominated by events viewed from off-axis directions, and the number of the 
off-axis events is about one-third of the on-axis emissions. We also compute the 
observed event rates of the GRBs, the X-ray-rich GRBs, and the X-ray flashes 
detected by HETE-2, and we find that they are similar. 

Subject headings: gamma rays: bursts — gamma rays: theory 



1. Introduction 

There is a strong correlation between the rest-frame spectral peak energy (1 + z)Ep and 
the isotropic equivalent 7-ray energy Eiso of the gamma-ray bursts (GRBs). This relation 
{Ep-E[sQ relation) was first discovered by Amati et al. (2002), and recently extended down 
to lower energies (Atteia 2003; Lamb et al. 2003b; Sakamoto et al. 2003), so that i^iso ranges 
over 5 orders of magnitude. A similar relation, the i?p-luminosity relation is also found 
by Yonetoku et al. (2003), and both relations could become a new distance indicator. The 
geometrically corrected 7-ray energies E^ = (1 — cosA6')i?iso narrowly cluster around a 
standard energy E^ ~ 10^^ ergs (Bloom, Frail, & Kulkarni 2003a; Frail et al. 2001), so 
that the opening half-angle of the jet in the on-axis uniform jet model ranges 2.5 orders 
of magnitude. This means that if the low isotropic energy events correspond to the wide 
opening half-angle jet, the jet opening half-angle of the typical GRBs A6 becomes less than 
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1° (Lamb, Donaghy, &: Graziani 2003a). However such a small angle jet has difficulties in 
the standard afterglow models and observations, (see also Zhang & Meszaros 2003). 

The low-energy (low-£'p) part of the relation consists of X-ray flashes (XRFs), that 
were identified by BeppoSAX (Heise et al. 2001) and other satellites (Strohmayer et al. 1998; 
Gotthelf et al. 1996; Hamilton et al. 1996; Arefiev, Priedhorsky, & Borozdin 2003) and have 
been accumulated by HETE-2 (Barraud ct al. 2003). Theoretical models of the XRF have 
been widely discussed (Yamazaki, loka, & Nakamura 2003b): high-redshift GRBs (Heise et 
al. 2001; Barraud et al. 2003), wide opening angle jets (Lamb, Donaghy, & Graziani 2003a), 
internal shocks with small contrast of high Lorentz factors (Mochkovitch et al. 2003; Daigne 
& Mochkovitch 2003), failed GRBs or dirty fireballs (Dermer et al. 1999; Huang, Dai, & 
Lu 2002; Dermer & Mitman 2003), photo sphere- dominated fireballs (Meszaros et al. 2002; 
Ramirez- Ruiz & Lloyd- Ronning 2002; Drenkhahn & Spruit 2002), peripheral emissions from 
collapsar jets (Zhang, Woosley, & Heger 2003) and off-axis cannonballs (Dar & DeRiijula 
2003) . The issue is what is the main population among them. 

We have already proposed the off-axis jet model (Yamazaki, loka, & Nakamura 2002, 
2003b). The viewing angle is the key parameter to understanding the various properties of 
the GRBs and may cause various relations such as the luminosity- variability /lag relation, the 
£'p-luminosity relation and the luminosity-width relation (loka & Nakamura 2001; Salmonson 
& Galama 2002). When the jet is observed from off-axis, it looks like an XRF because of 
the weaker blue-shift than the GRB. 

There are some criticisms against our off-axis jet model. The original version of our 

model (Yamazaki, loka, & Nakamiu^a 2002) required the source redshift to be less than 
~ 0.4 to be bright enough for detection, conflicting with the observational implications (e.g., 
Heise 2002; Bloom et al. 2003b). Yamazaki, loka, & Nakamura (2003b) showed that higher 
redshifts {z > 1) are possible with narrowly collimated jets (< 0.03 rad), while such small 
jets have not yet been inferred by afterglow observations (Bloom, Frail, & Kulkarni 2003a; 
Panaitescu & Kumar 2002; Frail et al. 2001). The luminosity distance to the sources at 
2;~0.4isdi~2 Gpc, which is only a factor of 3 smaller than that at 2; ~ 1 (corresponding 
to ~ 7 Gpc). Thus small changes of parameters in our model allow us to extend the 
maximum redshift of the off-axis jets to 2; > 1 even for ~ 0.1. This will be explicitly 
shown in § 3. Therefore, off-axis events may represent a large portion of whole observed 
GRBs and XRFs since the solid angle to which the off-axis events are observed is large. 

In this Letter, taking into account the viewing angle effects, we derive the observed 
Ep-Eiso relation in our simple uniform jet model. This paper is organized as follows. In § 2 
we describe a simple off- axis jet model for the XRFs. Then, in § 3, it is shown that the 
off-axis emission from the cosmological sources can be observed, and the Ep-Eiso relation is 
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discussed in § 4. Section 5 is devoted to discussions. We also show that the observed event 
rates of GRBs and XRFs are reproduced in our model. Throughout the Letter, we adopt a 
ACDM cosmology with {fl^, fl^, h) = (0.3, 0.7, 0.7). 



2. Prompt emission model of GRBs 

We use a simple jet model of prompt emission of GRBs considered in our previous works 
(Yamazaki, loka, & Nakamura 2003b; Yamazaki, Yonetoku, & Nakamura 2003). We assume 
a uniform jet with a sharp edge, whose properties do not vary with angle. Note that the 
cosmological effect is included in these works (see also Yamazaki, loka, & Nakamura 2002, 
2003a; loka & Nakamura 2001). We adopt an instantaneous emission, at t = to and r = Tq, 
of an infinitesimally thin shell moving with the Lorentz factor 7. Then the observed flux of 
a single pulse at frequency u = Uz/ (1 + z) and time T is given by 

dl Ml- (3 cos 9{T))r ' 

where 1 — f3 cos 6{T) = (1 + z)^^{cf3/ro){T — Tq) and Aq determines the normalization of the 
emissivity. The detailed derivation of equation (1) and the definition of A0(T) are found 
in Yamazaki et al. (2003b). In order to have a spectral shape similar to the observed one 
(Band et al. 1993), we adopt the following form of the spectrum in the comoving frame, 

f(u') = / (^V^o)'+"" expi-u'/u',) for u'/u', < as - Pb . . 

^ ^ I {i^'/O'^^Hc^B - /5b)"--^- expiPs - aB) for u'/u', > aB - Pb , ^ ' 

where i^'q, Ob, and Pb are the break energy and the low- and high- energy photon index, 
respectively. Equations (1) and (2) are the basic equations to calculate the flux of a single 
pulse. The observed flux depends on nine parameters: 7, a^, Pb, A6', Aq'J^, tq/Pc'j'^, ■jUq, 
z, and 9y. 



3. The Maximum Distance of the Observable BeppoSAX-XRFs 

In this section, we calculate the observed peak flux and the photon index in the energy 
band 2-25 keV as a function of the viewing angle 7^^. The adopted parameters are A^^ = 0.1, 
as = -1, Pb = -2.5, 7//^ = 300 keV and ro/Pc-/^ = 10 s (Preece et al. 2000). We fix the 
amplitude Aqj^ so that the isotropic equivalent 7-ray energy Eiso = 47r(i|(l + z)^^S{20~ 
2000 keV) satisfies the condition 

^{A9fEi,, ^E^ (= const.) , (3) 
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when 6y — 0. In this section, we take the standard energy constant — 1.15 x 10^^ ergs 
(Bloom, Frail, & Kulkarni 2003a). Then we obtain 7^Ao = 2.6 x 10^ erg cm^^. The redshift 
is varied from z — 0.01 to 1.0. 

For our newly adopted parameters and the spectral function in equation (2), we use a 
revised version of Fig. 3 in Yamazaki, loka, Sz Nakamura (2002), which originally assumed a 
different functional form of /(z^') and used the old parameters = 0.5 x lO^-*^ erg (Frail et 
al. 2001) and f3B = ~3. Figure 1 shows the results. Although qualitative differences between 
old and new versions are small, large quantitative differences exist. Since we now take into 
account the cosmological effects that were entirely neglected in the previous version, the 
observed spectrum becomes softer at higher z. The dotted lines in Figure 1 connect the 
same values of ^9^ with different z. The observed XRFs take place up to z ~ 1 in contrast 
to our previous result oi z — 0.4 and have viewing angles < < 2A6. The reason 
for this difference comes from the increase of the jet energy, the different spectrum and the 
different high-energy photon index. It is interesting to note that the only known redshift 
for XRFs so far is z = 0.25, one of the nearest bursts ever detected (Sakamoto et al. 2003; 
Soderberg et al. 2003). 

We roughly estimate the event rate of the XRF detected by MVFC/BeppoSAX (Ya- 
mazaki, loka, & Nakamura 2002). From the above results, the jet emission with an open- 
ing half-angle A^^ is observed as the XRF (GRB) when the viewing angle is within A^^ < 
6y < 2A^ (0 ^ ^ A^). Therefore, the ratio of each solid angle is estimated as 
/xrf//grb ~ (2^ — 1^)/1^ = 3. Using this value, we obtain i?xRF ~ 1 x 10^ events yr~^ 
for the distance to the farthest XRF c/xrf = 6 Gpc (see equation (5) of Yamazaki, loka, & 
Nakamura 2002). 

The derived value is comparable to the observation or might be an overestimation that 
may be reduced since the ffux from the source located at ^xrf ~ 6 Gpc is too low to be 
observed if the viewing angle 9^ is as large as ~ 2A6'. The ratio of the event rates of 
GRBs, X-ray rich GRBs (XR-GRBs) and XRFs detected by HETE-2 will be discussed in 
the following sections. 

4. E^-Ei^^ Relation of HETE Bursts 

In this section, we perform Monte Carlo simulations in order to show that our off-axis 
jet model can derive the observed Ep-Eiso relation and the event rate of the XRFs, the 
XR-GRBs and the GRBs detected by HETE-2. We randomly generate 10^ bursts, each 
of which has the observed flux given by equations (1) and (2). In order to calculate the 
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observed spectrum and fluence from each burst, we need eight parameters: ^ , as , , A^, 
^o7^('^o//9c7^)^ , ^u'q , z, and 9^. They are determined in the following procedure. 

1. We fix 7 = 100. The parameters a^, Pb, and A^^ are allowed to have the following 
distributions. The distribution of the low-energy (high-energy) photon index aB {13b) 
is assumed to be normal with an average of —1 (—2.3) and a standard deviation of 0.3 
(0.3) (Preece et al. 2000). The distribution of the opening half-angle of the jet, A^, is 
fairly unknown. Here we assume a power-law form given as /as (i(A^) oc (A6')~'' d{^6) 
for A^min < A^ < A^max- Wc take = 2 for the fiducial case, and adopt A^max = 0.3 
and A^min = 0.03 rad, which correspond to the maximum and minimum values inferred 
from observations, respectively (Frail et al. 2001; Panaitescu & Kumar 2002; Bloom, 
Frail, & Kulkarni 2003a). 

2. Second, we choose £^-y|^^Q°, which is the geometrically-corrected 7-ray energy of the 
source in the case of 2; = and 9^ — 0, according to the narrow log-normal distribution 
with an average and a standard deviation of 51 -|-log(1.15) and 0.3, respectively, for 
log (£'^ I ^^"5° /I erg) (Bloom, Frail, & Kulkarni 2003a). Then, the isotropic equivalent 
7-ray energy for z = and = is calculated as E-isofz=o — 2(A^)~^£^-y|^!^"^° to 
determine the flux normalization 74o7^(ro//3c7^)^. 

3. Third, we assume the intrinsic E^-Esso relation for 2; = and 9y — 0: 



This may be a consequence of the standard synchrotron shock model (Zhang & Meszaros 
2002b; loka & Nakamura 2002), but we do not discuss the origin of this intrinsic re- 
lation in this Letter. The coefficient ^ is assumed to obey the log-normal distribution 
(loka & Nakamura 2002), where an average and a standard deviation of log,^ arc set 
to —0.7 and 0.15, respectively We determine 72/9 such that the calculated spectrum 
vS^ has a peak energy £'p|^!j^° when ^„ = and z — Q. 

4. Finally, wc choose the source redshift z and the viewing angle 9^ to calculate the 
observed spectrum and fluence, and flnd E.^ and i^iso- The source redshift distribution 
is assumed to trace the cosmic star formation rate, and the probability distribution 
of 9v is P{9y) d9v oc sin6'^rf6'^. To determine the redshift distribution, we assume the 
model SF2 of the star formation rate given by Porciani & Madau (2001). 

We place a fluence truncation of 5 x 10"^ erg cm~^ to reflect the limiting sensitivity 
of detectors on HETE-2. Although the detection conditions of instruments vary with many 
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factors of each event (Band 2003), we consider a very simple criterion here. This fluence 
truncation condition is also adopted in Zhang & Meszaros (2003). 

Figure 2 shows a result. Among 10^ simulated events, 288 events are detected by 
HETE-2. The others cannot be observed because their viewing angles are so large that 
the relativistic beaming effect reduces their observed flux below the limiting sensitivity. 
Pluses (+) and crosses (x) represent bursts detected by HETE-2; the former corresponds 
to on- axis events {9.,, < A6) while the latter correspond to off- axis events (ft^ > A^^). The 
events denoted by dots are not detected. The numbers of on-axis and off-axis events arc 209 
and 79, respectively. Nearby events (-2^1) with large viewing angles can be seen. Such 
bursts are mainly soft events with {1 + z)Ep less than 60 keV. 

When 6y < A6, E^ is related to Eiso as E^ oc E^^ [see equation (4)]. The dispersion 
of pluses (+) in the Ep-Eis^ plane comes mainly from those of "intrinsic" quantities such as 
eXCo\ AO and {. 

1 /2 

On the other hand, even when 6^ > A9, the relation Ep oc E^J^ is nearly satisfied 

for the observed sources. The reason is as follows. For a certain source, as the viewing 

angle increases, the relativistic beaming and Doppler effects reduce the observed fiuence 

and peak energy, respectively. When the point source approximation is appropriate for 

the large-^^ case, the isotropic energy and the peak energy depend on the Doppler factor 

S = [7(1 - pcos{0^ - Ae))]-^ as Ei^^ oc S{20 - 2000 keV) oc and Ep oc 5, respectively 

(loka & Nakamura 2001; Yamazaki, loka, & Nakamura 2002; Dar & DeRujula 2003). Hence 

we obtain Eiso oc Ep~''°'\ Here {a) is the mean photon index in the 20-2000 keV band, 

1/2 

which ranges between 13b and as- Therefore we can explain the relation Ep oc E^J^ for 

(a) ~ aB ~ —1. On the other hand, when 9^ is large enough for Ep to be smaller than 

20 keV, we find i^iso oc E^~^b ~ E'^-^ or Ep oc E^^^ since (a) ~ I3b- In this case, the relation 

1/2 

deviates from the line Ep oc E.^^^ , and the dispersion of E^^o becomes large for small Ep. 



5. Discussion 

We have shown that our simple jet model does not contradict the observed Ep-Esso 
relation and extends it to lower Ep or Ei^^ values. The low-isotropic energy part of the 
relation is dominated by off-axis events. The number of off-axis events is about one-third 
of on-axis emissions. An important prediction of our model has been also derived, i.e., we 

1 /2 

will see the deviation from the present relation Ep oc E-^^^ if the statistics of the low-energy 
bursts increase. 
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HETE team gives definitions of the XRF and the XR-GRB in terms of the hardness ratio: 
XRFs and XR-GRBs are events for which \og[Sx{2 - 30 keV)/5^(30 - 400 keV)] > 0.0 and 
—0.5, respectively (Lamb et al. 2003b; Sakamoto et al. 2003). We calculate the hardness ratio 
for simulated bursts surviving the fiuence truncation condition, and classify them into GRBs, 
XR-GRBs, and XRFs. It is then found that all XRFs have redshift smaller than 5. The ratio 
of the observed event rate becomes Rgrb '■ Rxr-grb '■ Rxrf ~ 2 : 6 : 1. This ratio mainly 
depends on the value of q. When q becomes small, jets with large increases, and hence 
intrinsically dim bursts (i.e., low-E^isol^L'o^ bursts) are enhanced. Owing to equation (4), soft 
events are enhanced. In the case oiq—1 with the other parameters remaining fiducial values, 
the ratio is Rgrb '■ Rxr-grb '■ Rxrf ~ 1 : 9 : 3. For any cases we have done, the number 
of XR-GRBs is larger than those of GRBs and XRFs, while the event rate is essentially 
comparable with each other. HETE-2 observation shows Rgrb '■ Rxr-grb '■ Rxrf ~ 1 : 1 : 1 
(Lamb et al. 2003b). Although possible instrumental biases may change the observed ratio 
(Suzuki, M. & Kawai, N., 2003, private communication), we need more studies in order to 
bridge a small gap between the theoretical and the observational results. 

We briefly comment on how the results obtained in this Letter will depend on the Lorentz 
factor of the jet 7. If we fix 7 = 200, the relativistic beaming effect becomes stronger and 
less off-axis events are observed than in the case of 7 = 100; off-axis events are 13 % of the 
whole observed bursts when 7 = 200, while 27 % for 7 = 100. The ratio of the observed 
event rate for 7 = 200 is Rgrb '■ Rxr-grb '■ Rxrf ~ 2 : 5 : 1, which is similar to that for 
7 = 100. 

The Ep-Eiso diagram of the GRB population may be a counterpart of the Herzsprung- 
Russell diagram of the stellar evolution. The main-sequence stars cluster around a single 
curve which is a one-parameter family of the stellar mass. This suggests that the Ep-E^so 
relation of the GRB imphes the existence of a certain parameter that controls the GRB 
nature like the stellar mass. We have shown that the viewing angle is one main factor to 
explain the Ep-E^so relation kinematically. Our model predicts the deviation of this relation 
in the small £^iso region, which may be confirmed in future. 

In the uniform jet model, the afterglows of off-axis jets may resemble the orphan af- 
terglows that initially have a rising light curve (e.g., Yamazaki, loka, & Nakamura 2003a; 
Granot et al. 2002; Totani & Panaitescu 2002). The observed -R-band light curve of the 
afterglow of XRF 030723 may support our model (Fynbo et al. 2004). 
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0.1 1 10 100 

Fpeak^ 

Fig. 1. — Photon index in the energy range 2-25 keV as a function of the peak flux in the 
same energy range by varying the source redshift z. This figure is the updated version of 
Fig. 3 in Yamazaki, loka, & Nakamura (2002). We adopt 7A6I = 10, as = -1, /3b = -2.5, 
and = 300 keV. The values of the viewing angle 7^^ are given in parenthesis. Three solid 
curves correspond to 2; = 0.01, 0.1, and 1, respectively. The same values of 7^^, with different 
z are connected by dotted lines. The observed data of BeppoSAX-XKFs are shown from 
Heise et al. (2001). Squares (triangles) are those which were (were not) detected by BATSE. 
Two dashed lines represent observational bounds. Note that an operational definition of the 
XRF detected by Wide Field Cameras (WFCs) on BeppoSAX is a fast X-ray transient that 
is not triggered and not detected by the Gamma-Ray Burst Monitor (GRBM) (Heise et al. 
2001). In the region to the left of the vertical dashed hue, the peak flux in the X-ray band is 
smaller than the hmiting sensitivity of WFCs, and such events cannot be observed. In the 
region to the right of the oblique dashed line, the peak flux in the 7-ray band is larger than 
the hmiting sensitivity of the GRBM, and such events are observed as GRBs. 
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Fig. 2. — Distribution of simulated bursts in the (1 + z)Ep-Eiso plane. Pluses (+) and 
crosses (x) represent bursts that can be detected by HETE-2; the former is on-axis events 
{6y < AO) while the latter is the off-axis case {6^, > A6). The events denoted by dots are 
not detected. The dashed line represents the best fit of the observation given by E^ ~ 
95keV(Eiso/1052grgs)i/2 (Lamb et al. 2003b). 



